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ABSTRACT 

Aims. We investigate the nature of transverse kink oscillations of loops expanding 
through the solar corona and how can oscillations be used to diagnose the plasma pa- 



P4. 

_ rameters and the magnetic field. In particular, we aim to analyse how the temporal 

dependence of the loop length (here modelling the expansion) will affect the P1/P2 
£L| period ratio of transverse loop oscillations. 

Methods. Due to the uncertainty of the loop's shape through its expansion, we discuss 

5— 1 ' 

■ separately the case of the loop that maintains its initial semi-circular shape and the 
CZ2 . 

case of the loop that from a semi-circular shape evolve into an elliptical shape loop. 
The equations that describe the oscillations in expanding flux tube are complicated 
due to the spatial and temporal dependence of coefficients. Using the WKB approx- 
imation we find approximative values for periods and their evolution, as well as the 
OV period ratio. For small values of time (near the start of the expansion) we can employ 

\Q ' a regular perturbation method to find approximative relations for eigenfunctions and 

■ eigenfrequencies. 
Results. Using simple analytical and numerical methods we show that the period of 
oscillations are affected by the rising of the coronal loop. The change in the period due 
to the increase in the loop's length is more pronounced for those loops that expand 
into a more structured (or cooler corona). The deviation of periods will have signifi- 

J_l 1 cant implications in determining the degree of stratification in the solar corona. The 

effect of expansion on the periods of oscillations is considerable only in the process of 
expansion of the loop but not when it reached its final stage. 

Conclusions. The present study improves our understanding of the complexity of dy- 
namical processes in the solar corona, in particular the changes of periods of kink 
oscillations due to temporal changes in the characteristics of the coronal loop. Our 
results clearly show that the problem of expansion of coronal loops can introduce sig- 
nificant changes in the period of oscillations, with consequences on the seismological 
diagnostics of the plasma and magnetic field. 

Key words. Magnetohydrodynamics (MHD) — Sun: corona — Sun: magnetic fields — Sun: 
oscillations 
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1. Introduction 

Dynamical processes and transients observed in solar and space plasmas received consid- 
erable attention due to their ability to help scientists to diagnose remotely the parameters 
of the medium (temperature, scale-height, density, transport coefficients, etc.) in which 
these events occur, the magnitude and structure of the magnetic field, and the stability 
of the plasma. Seismological techniques and methods, imported from Earth's seismology 
and helioseismology, assume the combination of high resolution observations (amplitude, 
wavelength, propagation speed, damping time/length) with theoretical models (dispersion 
and evolutionary equation) in order to derive quantities that cannot be measured directly 
or indirectly. In particular, coronal seismology emerged as one of the most dynamically 
developing methods of solar physics (see, e.g. Roberts et al. 1984, Nakariakov et al, 1999. 
Ruderman and Roberts 2002, Andries et al. 2005, 2009, Ballai et al. 2005, Banerjee et 
al. 2007, Verth et al. 2007, Ballai 2007, Ruderman et al. 2008, Morton and Erdelyi 2009, 
Ruderman and Erdelyi 2009, Wang et al. 2012, just to name a few). 

One important candidate in coronal seismology are transverse kink oscillations, i.e. os- 
cillations which exhibit periodic movement about the loop's symmetry axis. Recent CoMP 
observations (Tomczyk et al. 2007) showed that the predominant motion of coronal loops 
is the transverse kink oscillation and this is the easiest to generate. The triggerring of these 
oscillations can be through a lateral forcing process at arbitrary height of the loop by a 
blast wave (or EIT wave) that emanates as a result of a sudden energy release by a flare 
and/or CME (see, e.g. Ballai 2007, Ballai et al. 2008). Secondly, kink oscillations can also 
be triggered by the transverse motion of the footpoints due to the granular buffeting of flux 
tubes in the photosphere. This scenario is true not only for coronal structures, but applies 
to all magnetic entities in the solar atmosphere that can serve as waveguides. Global EIT 
waves can also interact with prominence fibrils, as observed by, e.g. Ramsey and Smith 
(1966), and more recently by Eto et al. (2002), Jing et al. (2003), Okamoto et al. (2004), 
Isobe and Tripathi (2007), Pinter et al. (2008) (for a comprehensive review see, e.g. Arregui 
et al. 2012) in order to generate kink waves and oscillations in prominences. 

The dispersion relations for many simple (and some quite complicated) plasma waves 
under the assumptions of ideal magnctohydrodynamics (MHD) are well known; they were 
derived long before accurate EUV observations were available (see, e.g. Edwin and Roberts 
1983, Roberts et al. 1984) using simplified models within the framework of ideal and linear 
MHD. Although the realistic interpretation of many observations are made difficult by the 
spatial and temporal resolution of present satellites not being quite sufficient, considerable 
amount of information about the state of the plasma, and the structure and magnitude of 
the coronal magnetic field, can still be obtained. 

The mathematical description of waves and oscillations in solar structures is, in general, 
given by equations whose coefficients vary in space and time. It has been recognised by, e.g. 
Andries et al. (2005) that the longitudinal stratification (i.e. along the longitudinal symme- 
try axis of the tube that coincides with the direction of the magnetic field) is modifying the 
periods of oscillations of coronal loops. Accordingly, in the case of kink waves, these authors 
showed that the ratio P\jPi (where P\ refers to the period of the fundamental transverse 
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oscillation, while Pi describes the period of the first overtone of the same oscillation) can 
differ - sometimes considerably - from the canonical value of 2, that would be recovered if 
the loops were homogeneous. These authors also showed that the deviation of Pi /Pi from 
2 is proportional to the degree of stratification. This problem was also discussed in other 
studies such as Dymova and Rudcrman 2006, Diaz et al. 2007, McEwan et al. 2008, Ballai 
et al. 2011, Orza et al. 2012. In a recent analysis, Ballai et al. (2011) discusses about the 
ambiguity of the period ratio seismology, as some other effects could result in the observa- 
tion of multiple periods and each interpretation results in different value for the magnetic 
field and/or degree of stratification. 

Later, studies by, e.g. Verth et al. (2007), showed that it is not only density stratification 
that is able to modify the P1/P2 period ratio, but the variation of the loop's cross section 
area has also an effect on the period ratio. While the density stratification tends to decrease 
the period ratio, a modification of the cross section (i.e. when the magnetic field is flaring 
up as we approach the apex) tends to increase the Pi / Pi value. 

The investigation of properties of oscillations of coronal loops when equilibrium param- 
eters of the plasma depend on time is a relatively new area of dynamical studies in the solar 
atmosphere, in particular for coronal seismology. Morton and Erdclyi (209) and Morton et 
al. (2010) studied the effect of cooling (i.e. temporal dependence of temperature) on the 
dynamics of kink oscillations and travelling waves and they found that the cooling of the 
plasma results period decrease and amplification of oscillations and a mean to dissipate 
the energy stored in waves propagating in an unbounded plasma. The same idea was used 
later by Morton et al. (2011) to study the properties of torsional Alfven waves in coronal 
loops. The problem was recently re-considered by Ruderman (2011a) who showed that the 
cooling also generates an amplification of kink oscillations. The amplification of oscillations 
appears to be a competing effect with the damping due resonant absorption as shown by 
Ruderman (2011b) 

The problem of loop emergence and expansion through the solar atmosphere is one of 
the most challenging topics of solar physics as it involves the analysis of the evolution of the 
magnetic field in different regions of the solar interior and atmosphere where conditions can 
change from region to region. According to the standard theory, the magnetic field produced 
by the dynamo action in the tachocline is transported through the solar convective zone 
towards the solar surface by magnetic buoyancy coupled with convective motion (Parker 
1955, 1988). Once at the surface, the emerged flux tube creates sunspots and bipolar active 
regions (Zwaan 1987). In the solar atmosphere the rise of the flux tube continues due to 
an excess of the magnetic pressure inside the loop (e.g. Archontis et al. 2004). For the 
purpose of our investigation, we will assume that this excess is balanced at the transition 
region (TR) and from this height the expansion is not a driven problem any longer, instead 
the loop moves through the corona in the virtue of its inertia. During the emergence 
and expansion phase, the flux tube can interact with existing magnetic structure in the 
solar atmosphere and this might be responsible for the appearance of small-scale (e.g. 
compact flares, plasmoids, X-point brightcnings) and large-scale events (flares and CMEs) 
as suggested by Archontis (2004). Loop emergence often is associated with strong upflows 
as observed by, e.g. Harra et al. (2010, 2012). 
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It is rather straightforward to imagine what happens with oscillations in a loop when the 
length of the loop is increasing. As the length of the loop becomes larger, the frequency of 
oscillations becomes smaller, i.e. the periods of oscillations are expected to grow, however in 
an inhomogeneous waveguide particular periods will be differently affected by the combined 
effect of inhomogeneity and loop length increase. Therefore, we expect that the period ratio 
of oscillations will change not only with the degree of inhomogeneity but also with time. 

The aim of this paper is to investigate the effect of the loop expansion through the solar 
corona on the period ratio Pi / Pi and the consequences of the inclusion of the length of 
the loop as a dynamical parameter on estimations of the degree of density stratification. 
The paper is structured in the following way: in Section 2 we introduce the mathematical 
formalism and obtain analytical results for an expanding loop that initially starts from a 
semi-circular shape and this shape is preserved throughout the expansion. Later, in Section 
3 we generalise our findings by assuming that the expansion of the loop in the solar corona 
does not occur with the same speed in the horizontal and vertical direction, in this case, 
the initial semi-circular loop transforms into a loop with elliptic shape. Finally, our results 
are summarised in the last section. 

2. The mathematical formulation of the problem 

In our analysis we will capture the dynamical behaviour of the loop once reached the 
transition region (TR). At this height, the loop can interact with a blast wave propagating 
in the low solar corona. This instigator could be easily identified with an EIT wave that 
propagates in the low corona over very large distances and they are known to be one of 
the major source for kink oscillations of coronal loops (see, e.g. Ballai 2007). We assume 
that the height at which the loop starts its journey through the corona is at 3 Mm above 
the solar surface. A typical loop length is of the order of 300 Mm. For practical reasons, 
we are going to consider that the height of the loop in its final position would be about 97 
Mm resulting in a loop length of about 305 Mm. 

The raising speed of loops is generally taken to be 10-15 km s _1 (Chou and Zirin 1988, 
Archiontis 2008), therefore the time needed for the loop to travel the distance from the 
TR to its steady position is easily estimated to be between 3.4 and 5.2 hours. This time is 
at least two order of magnitude larger than a typical period of kink oscillations, so there 
is enough time for the development of oscillations. Higher rising speeds are also possible, 
recently Schmidt and Ofman (2011) reported expansions of a post-flare loop with speeds 
of hundreds km s _1 . Standing waves are formed if the speed of change in the length of 
the loop is smaller than the period of oscillations. This condition is easily satisfied for fast 
kink oscillations. We assume that the loop expands into the "empty" corona, i.e. it will not 
encounter any interaction with existing magnetic elements. 

In the first instance we assume that the semi-circular shape of the loop at the TR is 
preserved throughout the expansion. Assuming an initial height of 3 Mm above the surface, 
the distance between the footpoints of the loop is 6 Mm. When reaching the final height 
of 190 Mm, the footpoints travel over a distance of 197 Mm. In addition we assume that 
the expansion occurs at constant temperature (isothermal process) and throughout the 



4 



Ballai and Orza: Transverse kink oscillations of expanding coronal loops 

expansion the loop maintains its cross-section constant. Due to the increase in the volume 
of the loop, a pressure difference is generated meaning that the plasma flows along magnetic 
field lines resulting in a density that depends not only on the height, z, but also on time. As 
it raises through the solar corona, the tangent to the loop is also changing monotonically, 
therefore the momentum equation in an equilibrium state becomes 
du du dp 

+ Pqu— = — gpocos/3, (1) 

at oz oz 

where the angle /3 is a function depending on time and space, u is the equilibrium flow of 
the plasma, po & n d Po are the equilibrium density and pressure, and g is the gravitational 
acceleration. In coronal loops, the flows are of the order of a few tens of km s _1 , therefore 
in Eq. ((1} we can neglect the terms on the left hand side since 

du u 2 10 8 

u d-z~T~W = 0{1) > 

that is two order of magnitude smaller than terms on the right hand side. We can write 
pressure as 

, ,n k B T p (z,t) 
Po{z,t) = 

TO 

where k B is the Boltzman constant, To is the constant temperature and to is the mean 
atomic mass per particle. Introducing this expression into the RHS of Eq. ([!} we obtain 
that 

1 dp Q cos /3(z,t) _ 

Po dz H ~ [2) 

where H is the constant density scale-height. Focussing on the density distribution inside 
the loop, we can integrate the above equation to obtain 



pj exp 



cos p{z',t) ; 

H 



(3) 



Assuming that the loop is semicircular we obtain that /3(z,t) = nz/L(t), therefore the 
density inside the coronal loop becomes 



pj exp 



Lit) 7TZ 

■ sm ■ 



(4) 



nH L(t)_ 

where pf is the density of the plasma at the footpoint of the loop. For simplicity we will 
assume that the external density can simply be written as p e = Dpi, where throughout 
our calculations we consider D — 0.5. Figure 1 depicts schematically the change of the 
equilibrium density, both in space and time. Here length was normalised to the length of 
the loop at the start of the emergence (considered Lq — 3tt Mm) and time (here denoted 
by r) was normalised to the quantity Lo/v r , where v r is the vertical rising speed of the 
loop considered v r = 15 km s _1 . According to our expectations, the density of the loop 
decreases with time. We assume that the flux tube is thin even at the beginning of its 
expansion. 

In the thin flux tube approximation the dynamics of transverse kink oscillations is given 
by (see, e.g. Ruderman 2010, 2011a) 

d TT d V ( d TT d V 2B 2 d 2 r 1 



Jt + Ui 8-z) ^U + C ^J r? -^ = °' (5) 
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Density profile 




Fig. 1. A schematic representation of the evolution of the equilibrium density measured on 
the vertical axis in the units of density at the footpoint. Here lengths are given in units of 
the loop length at the start of the expansion in the corona (Lo) and time is given in units 
of L /v r , where v r is the constant rising speed in the vertical direction, here taken to be 
15 km s _1 



where Ui and U e are the internal and external flows, B is the magnetic field strength 
(here assumed to have identical value inside and outside the loop), p is the magnetic 
permeability of the free space and 77 is a complex valued displacement of the loop, with 
77 = tjji + ir] 1. In Cartesian coordinates the loop's displacement in the x and y directions, 
£ x and £ y , are given by 

£x = VR, = ~Vi- 

We assume a quasi-stationary equilibrium, so that the characteristic time variation of 
equilibrium quantities (t c h) is much longer than the period of kink oscillations (P) and 
introduce the small parameter e so that P = et c h- Using the definition of the period of 
oscillations we can write 

L(iJoPch) 1/2 , ^ R -i L(nop ch )^ 2 

= etch B = e , (b) 

-t>0 tch 
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meaning that we can introduce a scaled magnetic field, so that Bo = (Bq. As a result, the 
equation describing the dynamics of the kink oscillations can be written as 

d TT d\ 2 ( d TT d\ 2 2e- 2 B 2 d 2 V n 

The above equation must be solved subject to the standard boundary conditions, i.e. 
17(2: = 0,z = L) = 0. 

Following the solution method proposed by Ruderman (2011a) we will solve Eq. ((?]) using 
the Wentzel-Kramers-Brillouin (WKB) method (see, e.g. Bender and Orszag 1987) and 
assume that the solution of the equation will be of the form 



77 = ^~VS fc (z,t) exp 

fc=0 



e 



(8) 



In the first order of approximation (often called geometric optics), Eq. ([7]) reduces to 



where 

d$(t) „ 2 2B 2 



dt ' fl (pi(z,t) +p e (z,t))' 

Equation © must be solved subject to the boundary condition So — when z — and 
z = L. Equation Q together with the line-tying condition form an eigenvalue problem in 
which n is the eigenvalue and ft 2 is a real function. 

In the next order of approximation (also called the approximation of physical optics), 
Eq. ([7]) reduces to 



d 2 Si w_ s _ mi 

dz 2 d 2 K 1 ~ d 2 K 



dS S dQ piUi + p e U e dSo 



(10) 



dt 2Vt dt pi + p e dz 

that has to be solved subject to the boundary condition 6*1(2 = Q,z = L) = 0. The 
boundary-value problem determining S\ has a solution only when the RHS of Eq. (|10j) 
satisfies the compatibility condition, i.e. the orthogonality to So- After multiplying the 
RHS of Eq. (fTU]) by So and integrating with respect to z in the interval (0, L) we obtain 
(similar to Ruderman 2011a) that the compatibility condition reduces to 

f L S 2 

lo I —^-dz = const., (11) 

Jo C K 

where 

lj = e~ fl, ck — €~ 1 5k- 

As a consequence, the dynamics of kink oscillations in coronal loops is fully described by 
the system of equation ([9]) and (|TT|) . In deriving Eq. ((TTj) we took into account the mass 
conservation equation, relating the plasma flow and its density 

dp + d(pu) = Q 

dt dz 

Let us now discuss a special case when the density depends on time only. This case 
would correspond to an initial expansion of the loop when the height of the loop is less 
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than the scale-height (assuming expansion into an isothermal 1 MK corona, this height 
would correspond to 47 Mm), In this case, we may expect that the amplitude of oscillations 
increases as the loop expands. Indeed, it is easy to show that the amplitude of oscillations 
behaves like 

where A(0), Ck-(O) and Lq are the amplitude of oscillations, the kink speed and the length 
of the loop at t = 0. 

A particular case to discuss is when the expansion of the loop occurs linearly with time 
and we write that 



L(t) = L + v r t 



(13) 



where Lq is the length of the loop at the initial time, i.e. at the TR level (Lq = 3tt Mm) 
and v r is the rising speed, here assumed constant. Accordingly, the kink speed becomes 



4 



c Kf exp 



L + v r t 



Lq + V r t 



(14) 



where 



2B 2 



' Ki (i oP f(i + Dy 

Returning to the general case, the equations describing the dynamics of kink oscillations 



—S-dz = const. 



are 

d 2 S ^ 

OZ 1 CK 1 Jo 

Let us introduce a new set of dimensionless quantities 

v r t uLq 

= — , 1.1 : 

^0 



t-j- 



CKf 



n 

^0 — 

^0 



In the new variables the equations to be solved transform into 
d 2 S 



and 



exp 



-L (l + r) . < 
■ sin ■ 



ttH 



1 



Q Z S = 0, 



Sq exp 



-L (1+t) . tt^ 
■ sin ■ 



ttH 



1 



e?£ = const. 



(15) 



(16) 



(17) 



(18) 



that should be solved subject to the boundary conditions 5o(£ = 0; £ = 1 + r) =0. The 
solution of the system (fTTj) - ([T5]) can be found numerically using, e.g. the shooting method. 
In Fig. 2 we display first the variation of periods of oscillations of the fundamental mode 
and its first harmonic for three different values of H in terms of the dimensionless time 
variable, r. The bands for each period are clearly labelled in the figure. The three distinct 
value of periods were obtained for three values of scale-height keeping the initial length 
of the loop at 3ir Mm. The dotted line corresponds to a scale-height of 70.5 Mm, which, 
assuming a plasma in hydrostatic equilibrium, would correspond to a plasma temperature 
of 1.5 MK. The solid line is plotted for a loop expanding into a corona where the constant 
scale-height is 47 Mm, that would correspond to a 1MK hot plasma. Finally, the dashed 
line stands for an expansion of the loop into a plasma where the density scale- height is 23.5 
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Fig. 2. The variation of the periods of the fundamental mode and its first harmonic with 
the dimensionless time variable t for three different values of stratification: H=70.5 Mm 
(dotted line), H=47 Mm (solid line), and H=23.5 Mm (dashed line). 



Mm, corresponding to a plasma temperature of 0.5 MK. The two bands for the periods 
clearly show that the two oscillations are differently affected by the expansion, i.e. change 
in the length of the loop. This is also obvious in Fig. 3 where we plotted the ratio of 
the periods of the fundamental and first harmonic as a function of the time variable, r. 
The plot clearly shows that the oscillations of a loop expanding into a "hot" plasma (i.e. 
large scale-height) are the least affected, but in all three cases the period ratio decreases 
with time. We should note here that the periods shown in Fig 2 do not start at the value 
of 2 because at the start of their expansion through the solar corona, they are already 
stratified and the least stratified is the case that corresponds to H = 1 .5 MK. Let us make 
a final note: in observations, the identification of periods in coronal loops is a dynamical 
process, i.e. in EUV the intensity is measured in one location (or mega-pixel to reduce 
errors) for a long time-period. The duration of observation vary, and mostly they are 
driven either by the availability of the instrument or by the limited life-time of oscillations 
before they are damped by, e.g. resonant absorption, although observations show that not 
all loop oscillations damp (probably these oscillations are maintained by a constant lateral 
buffeting). Typically, for damped oscillations, the detection time is a few periods with a 
range of 6.7-90 minutes (Aschwanden 2004). In terms of the dimensionless quantity t, this 
range would correspond to r = 0.64 — 8.6. Although the lower limit is too small to count in 
the effect of expansion, a duration of r = 8.6 would add an important effect in studying the 
transverse kink oscillations. Fig. 3 shows that for one given value of internal structuring, 
the value of the period ratio can change also because of the expansion of the loop. Although 
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P1/P2 for t = (0.0,10.0) 



1.95 - 



1.90 - 



1.85 - 



1.80 - 



1.75 




Fig. 3. The variation of the P1/P2 period ratio with respect to the dimensionless parameter 
r, for a loop expanding in the solar corona with persisting semi-circular shape. The meaning 
of each line-style is identical to Fig. 2. 



the periods shown in Fig. 2 display a monotonic increase with the parameter r, in reality 
these values will saturate, the saturation occurring faster for those modes that propagate 
in a highly structured plasma (e.g. for the case of H=23.5 Mm, the saturation value of 
periods is about 400 seconds and the saturation starts at about r = 20). In addition, the 
period ratio for all cases discussed here tend to same value (near 1) for large values of r. We 
need to mention here the very important fact that the effect of expansion on the periods 
of oscillations and the period ratio is relevant only in the expansion phase of the loop. 

Analytical solutions of the Eqs. (ITTt can be obtained for the limiting case of a loop at the 
beginning of its expansion through the solar corona, i.e. small values of r (see Appendix). 
The results confirm the tendency of periods to increase with time and of the period ratio 
to decrease with time. 



3. Non-circular emergence 

In reality, the expansion of a loop in the empty corona so that the semi-circular shape is 
preserved is unlikely since the foot points have to move in a much denser plasma than the 
apex of the loop. That is why in this section we will assume that the expansion rate in 
the vertical direction is larger than the expansion of footpoints in the horizontal direction. 
The expansion still starts at the TR level where the shape of the loop is semi-circular and 
assume that the process remains isothermal. As a result of different expansion rates, the 
loop evolves so that the shape becomes more elliptical. The properties of transverse loop 
oscillations in an elliptical coronal loop was studied recently by Morton and Erdelyi (2009) 
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assuming that the semi-elliptical shape is reached in the emerging stage, before reaching 
a semi-circular shape. Although their topic is related to the research of the present study, 
the problem of expansion is a dynamical process that should be treated accordingly. They 
found that the difference in P± / P2 period ratio between the circular and elliptical shape is 
up to 6%. 

Since the dynamics is going to be different over the two directions, it is more convenient 
to introduce a polar coordinate system in which 

x = a(t) cos 6, z = b(t)sm9, (19) 

with a(t) < bit), the overdot denotes the derivative with respect to time and the length 
of the loop is covered by the parameter 9 that varies now between and tt. It is more 
convenient to use the coordinates along the loop, s, therefore the dynamics of transverse 
kink oscillations is described by 

8 8\ 2 (8 8\ 2 2B 2 8 2 n 

di + ^dl) ^ + ^{dt +U ^s) *~^a?= ' (20) 
A key parameter in our discussion is going to be the arc-length that is defined as 
8s 



^ = Ja(t) 2 sin 2 6 + bit) 2 cos 2 6 = a(0, t). (21) 
06 v 

In order to solve Eq. (|2"0|) we would need to express the density as function of s and t. 
However, it turns out that it is much easier to deal with the variable 9 instead. Therefore 
we express the derivatives in the governing equation as 
d__ld_ 1 d 2 sm29(a(t) 2 - b{t) 2 ) 8 

ds~ad6' 8s 2 ~a 2 88 2 ~ 2a 4 09' ( ' 

Assuming again a quasi-stationary equilibrium , similar as in the previous section, and 
introducing the small parameter, e, the governing equation for transverse kink oscillations 
can be written as 



8 2 V , 2 f p i U i + Pe U e \ d 2 rj , 1 



dt 2 a \ pi + p a I d9dt a 2 
d 2 rj 



Pi U 2 + Pe U 2 

Pi + Pe 



2e- 2 B 2 



Po(Pi + Pe) 



89 2 



PiUi + p e U e aa sin 9 + bbcos 2 9 

2 n 



| p t U 2 + Pe U 2 sm28(a(t) 2 - b(t) 2 ) 
P i + p e 2a 4 

er 2 B 2 sin28 (a(t) 2 - b(t) 2 ) 
Po(Pi+Pe) 2a 4 

that has to be solved subject to the boundary conditions rj(9 = 0, 9 = %) = 0. 

Again, we will solve this equation using the WKB approximation presented earlier and 
suppose a solution identical to Eq. ((HJ. In the first order of approximation we obtain the 



|f=0, (23, 



So = 0, (24) 



equation 

d 2 S Q _ sm28(a 2 -b 2 )dS a 2 tt 
89 2 ~ 2a 2 89 + al- 

together with the usual boundary conditions at the two ends of the loop. In the next order 
of approximation we obtain 

d 2 S x sm28(a 2 ~ b 2 ) dSi a 2 fl 2 
d9 2 ~ 2a 2 89 + c 2 K 1 ~ 
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2iaVt 



K 



dS aS dfl pjUj + p e U e dS 
dt ' 1 m~dt Pi + p e W 



(25) 



This equation has to be solved subject to the boundary condition S\{6 = 0, 6 = tt) = 0. 
Again, the equation for Si will have solution if the right-hand side of the above equation 
satisfies the compatibility condition, i.e. the orthogonality to Sq. Following the same solving 
procedure as presented earlier, it is easy to show that the condition reduces to 

-g-d9 = (26) 

C K 

Therefore, the system of equations (j2"4"]l and (|2"6"|) will determine completely the dynamics 
of the expanding coronal loop. 

Let us assume that at t = the loop is semi-circular and its length is Lq. In order 
to reproduce the different movement over the two directions, we introduce two different 
expansion speeds in the horizontal (vh) and vertical (v v ) direction, so that u/, < v v . Again, 
we suppose that the motion occurs linearly in time and write the dynamics over the two 
axes as 

x = ^— + v h tj cos 9, z = ^— + v v tj sm9. (27) 
Let us introduce a new set of dimensionless quantities 

r=— , u= , a=—, U=—, S = —. (28) 

Lq CKf L V v Lq 

In the new variables, the governing equations become 



d 2 S , _ 2 _ 2 [ L (l + T7r) 



dff 2 



a lo exp 



sm ( 



Sq- 



sin 9 cos 9[t 2 it 2 (1 - U 2 ) + 2ttt(1 - U)} 
1 + T 2 7r 2 (sin 2 6 + U 2 cos 2 6) + 2r7r(sin 2 9 + U cos 2 6) X 

§ = »• w 

and 

= 0, (30) 

C K 

where ck = ck/v v . Figure 4 displays the evolution of the period of oscillations for the 
fundamental mode and its first harmonic for three different values of Lq/H, similar values 
as used in the previous section. Comparing the findings in Figs 2 and 4, the effect of 
the expansion into an elliptical shape compared to the constant semi-circular shape is 
evident. As time progresses the period of oscillations tend to a higher value for elliptical 
shape, however this conclusion is more true for the fundamental mode. The period of the 
fundamental mode corresponding to an expansion in a solar corona where scale height 
is only 23.5 Mm (corresponding to a temperature of 0.5 MK in a loop in hydrostatic 
equilibrium) saturates rather quickly. 

A more significant change is evident when comparing the P1/P2 period ratio of the 
expansion into an elliptical shape shown in Fig. 5 to the variation of the period ratio for 
a semi-circular shape. For the same time interval the decrease of the period ratio is much 
more significant, and similar to the previous case, the period ratio is more affected for the 
case of strong stratification, i.e. small H. 
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P1 and P2 for t = (0. 0,10.0) 
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Fig. 4. The same as in Fig. 2, but here we assume that the expansion of the loop occurs 
such that the loop evolves into a loop with a semi-elliptical shape. The meaning of different 
line-styles is identical to Fig 2. 



P1/P2 for t = (0.0,10.0) 




Fig. 5. The same as in Fig. 3, but here we assume that the expansion of the loop occurs 
such that the loop evolves into a loop with a semi-elliptical shape. The meaning of different 
line-styles is identical to Fig 2. 
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4. Conclusions 

The solar corona is a very dynamical environment where changes in the dynamical state of 
the plasma and field occur on all sort of time scales. In the present study we combined for 
the first time two kinds of dynamical events: the time-evolution of a coronal loop through 
its expansion into the "empty" corona and the transverse kink oscillations of coronal loops. 
The emergence and expansion of a coronal loop is a very complex phenomenon, but here we 
reduced our model to a simplified process, where the expansion is solely described by the 
change in the length of the loop with an associated temporal equilibrium density variation. 

The governing equation for kink oscillations was solved in the WKB approximation 
when the boundary conditions are time-dependent. As expected, due to the change in the 
length of the loop, the amplitude and periods of oscillations increase with time, however, 
the period ratio of the fundamental mode and its first overtone decreases. This last physical 
parameter is of paramount importance for the remote determination of density structuring 
of coronal loops with the help of seismological approaches. In the first instance we regarded 
the loop to have an initial semi-circular shape that is maintained through the expansion 
phase. Later, this restriction was lifted based on the natural assumption that the expansion 
into the vertical direction (i.e. in the direction of density decrease) occurs much easier than 
in the horizontal direction. In this limit, the loop evolves into a semi-ellipse, with the major 
axis in the vertical direction. Comparing the results of the two approaches it is clear that 
the behaviour of the period ratio is rather sensitive to the geometrical shape of the loop, 
a more significant drop in the P1/P2 ratio being achieved in the second case. Although 
our numerical result were obtained for three different structuring degrees (measured by the 
ratio of the initial loop length to the density scale-height) it is also evident that both the 
temporal change in the loop length and the stratification will have the same effect upon 
the period ratio resulting in a mutual amplification of the effect. 

Our model predicts that the amplitude of oscillations increases with time, however due 
to the particular choice of density, damping processes were neglected. Once the density is 
allowed to vary also in radial direction, according to the theory of resonant absorption (e.g. 
Goossens et al. 1992, Rudeman and Roberts 2002), loops will damp very quickly with the 
resonant position displaying a steady motion due to the change of the length of the loop. 
The amplitude of oscillations can also be damped due to the cooling of the plasma (Morton 
ct al. 2010), an effect that was also neglected here. In an expanding loop, the growth of the 
amplitude due to emergence and decay of amplitude due to resonant damping or cooling 
will be competing processes and the competition between these two effects will be discussed 
in a forthcoming paper. 

Acknowledgements. IB acknowledges the financial support by NFS Hungary (OTKA, K83133). The au- 
thors are grateful to Prof M.S. Rudemman and Dr. V. Archontis for their help and comments. 
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Appendix A: Solutions to the wave equation in the case of the loop at the 
beginning of its expansion 

An interesting insight into the character of the solution of the governing equation can be obtained analyt- 
ically in the limiting case of 

C = L (1 + t)/ttH < 1, 

i.e. we restrict ourself to the first part of the emergence into the solar corona. In this case the argument of 
the exponential function in Eq. 1171 can be expanded (keeping only the first two terms) and the equation 
to be solved reduces to 

— — + ui S = - sin — — CSb, (A.l) 
1 + r 

Since we are looking for periodic solutions and we expect that both the amplitude and frequency will 
depend on time, we will employ the Poincare-Lindstedt method (Meirovitch 1970) to find corrections to 
the eigenfunctions and eigenfrequencies (it can be easily shown that this method is similar to the re- 
normalization technique used by Ballai et al 2007 in the case of dispersive shocks). We are looking for 
solutions in the form of series and write 

oo oo 

So = ^C fe ^ fc) , ^ = ^C fc - (fc) , (A.2) 
fc=l fc=l 

where the functions Si are periodic functions. We first concentrate on the fundamental mode. After inserting 
the expansions IIA.2II into the governing equation lA.ll l , we collect terms proportional to subsequent powers 
of In the first order of approximation, collecting terms ~ 0(£°), results in 



9 2 S, (0) 



^-+^('») 2 S f ( ) 0) =0. (A.3) 



<9£ 2 

Solving this equation subject to the mentioned boundary conditions, yields in the case of the fundamental 
mode 

5(°>~ sin ^-, a>t°) ^— . (A.4) 

1+T 1 + T ^ ' 

In the next order of approximation (i.e. terms ~ 0(f)) we obtain 
. inn ~M "i 7r 2 . t 7r£ 



-2- + Z^S^ = sin 2 

de 2 (l + r) 2 1 + T 

2n ~(1) 



iVWsin — — (A.5) 

1+T 1+T 

that has to be solved subject to the boundary conditions 

^ 1) « = 0,5 = l + »?)=0. (A.6) 

The last term in Eq. I IA.5I I will cause secular growth of the solution, therefore rendering the solution Si 
non-periodic. To suppress this possibility we choose Cb\ = 0. As a result, the solution of Eq. IIA.5I) together 
with the boundary conditions is simply written as 

6A ~ sin — cos 3 + cos (A. 7) 

l + r 2 1 + T 6 V 1 + T/ v ' 

In the next order of approximation we collect terms C>(f 2 ) and obtain 
^ + ^ S f = -*<°) 2 S« sin ( J*-) - 

a? \i + tJ 

2^°>d)( 1 >S 6 )) . (A.8) 

Using the expression of u>(°\ Sq ' and Sq 1 ^ determined earlier, the RHS of the above equation can be 
written as 

t 2 2W 2 )~ 



RHS ■ 



3(1 + t) 2 1 + t 



sm h A sm 

l + T 1 + 7 
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2tt£ 3tt£ 
B sin h C sin ■ 



1 



1 - 



(A.9) 



where the coefficients of higher harmonics (A, B and C) are not needed for our discussion. In order to 
prevent non-periodic behaviour we need to impose the condition that the coefficient of the first term is 
zero, leading to 



5(2) 



6(1 +t) 

As a result, the eigenfunction and eigcnfrcqucncy of fundamental mode oscillations can be written as 



(A.10) 



5 , f 

bo = sin h ( 



3 tt£ 
— cos 

2 1 + 7 



2tt£ 
1 + r 



(A.ll) 



6(1 + r)' 

meaning that the change in the frequency due to expansion is a second order effect. It is easy to show that 
the period of the fundamental mode in this approximation can be written as 
12(1 + t) 



Pi 



(A.12) 



6 + C 2 

proving the increase of the period P\ with r seen in Fig 2. 

Repeating the same method for the first harmonic, where the eigenfunction and eigenfrequency in the 
zcroth-order approximation are 
27t£ _, m 2tt 



Sq = sin 



we obtain that 
2tt 



+ C^ 



-,(o) 



2tt 



(A.13) 



1 + r ■ 3 15(1 + r)' 
meaning that the period of the first harmonic in this approximation behaves like 

showing an increasing tendency with respect to r. Now, using Eqs. (I A .121) and l lA.1411 we can calculate the 
period ratio of the fundamental mode and its first harmonic as 



Pi 
Pi 



2 1 



60 



L 2 (l + r) 2 
60tt 2 H 2 



(A.15) 



so that the change in the period ratio is very small but decreases with f . 
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